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SUMMARY 

Acoustic measurements in air-water mixtures at a fixed temperature 
and pressure in which void fraction is varied yield a 
sound speed curve as a function of void fraction exhibiting 
a minimum in the vicinity of a void fraction of one half. A metastable 
theory in which the entropies of the mixture , gas phase , and liquid phase 
are individually held constant yields a corresponding curve whose mini- 
mum sound speed is about 20$ higher than the measured value. When an 
ad hoc modification is made which replaces isentropic by isothermal 
propagation the new theoretical minimum region is brought into agree- 
ment with experiment; however, the boundary values of the sound speed 
at void fractions of zero and unity are no longer the usual isentropic 
sound speeds for pure water and pure air. In this paper a new metastable 
theory is developed. By imposing only constant mixture entropy and phase compo- 
sition during propagation, the new theory successfully predicts both the minimum 
region and the end points. The theory produces a heat capacity ratio for the 
mixture which self-adjusts to both normal values for the pure phases and 
to values near unity for a broad range of void fraction values about the 
minimum in the curve. The new theory is extended to single-component pure 
substances including para-hydrogen and nitrogen with emphasis on the latter. 


^Associate Professor of Mechanical Engineering, Aeronautical Engineering, 
and Mechanics, Rensselaer Polytechnic Institute, Troy, New York. 



INTRODUCTION 


A large transonic wind tunnel is under construction at the Langley 
Research Center. This new National Transonic Facility (NTF) will use a 
relatively novel concept to achieve high Reynolds numbers - liquid nitrogen 
is injected directly into the nitrogen test gas and the resulting evaporation 
of the liquid cools the test gas to cryogenic temperatures approaching 
liquid-vapor saturation temperatures (see references 1-3). In support of 
this new wind tunnel, studies have been undertaken at Langley to insure 
that the low temperature behavior of the nitrogen test gas is sufficiently 
understood to guarantee its usefulness for transonic aerodynamic testing 
(see references k and 5). 

The present study involves -the calculation of sound speed for a gas-liquid 
mixture. In cryogenic tunnels such as the NTF, two-phase flow is always 
present at the liquid nitrogen injection station, and at times may be desired 
in the test section in the form of fog droplets for the purpose of flow 
visualization or for seeding for laser velocimetry. While the gas-liquid 
mixtures in cryogenic nitrogen tunnels are expected to involve mostly 
gas - vapor mass fraction generally greater than O .98 at the injection 
station and close to 1.00 in areas of seeding - the general sound speed 
literature used in this paper deals with mixtures of air-water and steam- 
water that are mostly liquid. Nevertheless, the literature can be used 
to develop the theory that will apply to cryogenic nitrogen tunnels. The 
present paper will consider these water mixtures before addressing . 
gas-liquid nitrogen. All nitrogen calculations will use thermodynamic 
properties from reference 6. 
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equilibrium, zero-frequency, two-phase, sound speed, m/sec 
metastable, two-phase sound speed, defined by equation (32), m/sec 
metastable, ad-hoc, two-phase sound speed, defined by 
equation ( 7 ), m/sec 

metastable, triply isentropic, two-phase sound speed, 
defined by equation (4), m/sec 

constant volume specific heat capacity, J/(kg-K) 

constant pressure specific heat capacity, J/(kg-K) 

2 

pressure, N/m 

specific entropy, J/(kg-K) 

non-dimensional correlation parameter defined in equation (39) 
specific gas constant, J/(kg-K) 
temperature, K 

3 

specific volume, nr/kg 

vapor mass fraction (quality) 

value of quality for which (3s/3T) is zero 

x 

void fraction j fraction of vapor volume to the total volume 
specific heat ratio (Cp/C^.) 

O 

density, kg/m 


Subscripts 

c pertaining to the critical point 

G denoting gas (vapor) phase 
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denoting liquid phase 
denoting isothermal value 



BACKGROUND 


Propagation of sound in two -.phase, gas-liquid systems has under- 
gone extensive theoretical and experimental study. One of the earliest 
expressions for this acoustic velocity comes from the so-called isen- 
tropic, homogeneous equilibrium model (IHE) in which average equilibrium 
properties are used in conjunction with the assumption of total phase 
equilibrium to describe sound propagation when the two-phases are finely 
interspersed in one another. In this paper, this velocity of sound is 
denoted by a Q and it is computed only for nitrogen. Details of the 
a Q computation are omitted because, first, it is such a widely known 
model and, second, the subsequent main thrust of the paper will deal 
with metastable models. As usual, the square of a^ is taken as (9P/9p) g 
and this quantity is evaluated herein at equilibrium two-phase states 
for nitrogen using a computer program generated at the Langley Research 
Center from reference 6. 

One of the chief characteristics of a is associated with the 

o 

manner in which an isentrope crosses the phase boundary as it passes 
from the one-phase region (liquid or vapor) into the two-phase region of 
a P, P property diagram where the saturated liquid and vapor curves 
are distinct and separate (in contrast to a P, T plane of pro- 
perties where the saturated liquid and vapor curves coalesce into a single 
curve). In the P, p plane the isentrope crosses a phase boundary with 

discontinuous slope; that is (9P/9p) has two values at a point on the 

s 

phase boundary. If this partial derivative is interpreted as the sound 
speed, then, from a total equilibrium point of view, there are two values of the 
sound speed at a point of the phase boundary. (As will be seen 
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later, it is possible to avoid having two sound speeds at a single point 
of the phase boundary by introducing metastable effects into the two- 
phase propagation of sound. The two-phase, equilibrium (3P/8p) is 
then simply a thermodynamic partial derivative to which no further physi- 
cal attribute is assigned.) Furthermore, the formula for a Q will not 
yield the usual single-phase velocity of sound at either the liquid or 
the vapor boundary. If the boundary values of a are denoted by a 

O U j Jj 

on the liquid side and a r on the vapor side and if the usual single- 
phase sound speeds are denoted by a^ and a^ - all in meters per 
second - then this discontinuous boundary behavior for the IHE model 
can be illustrated for nitrogen in the table that follows : 

T SATURATED LIQUID LIMIT SATURATED VAPOR LIMIT 


(K) 

ONE-PHASE 

a L 

TWO-PHASE 

a T 
O , L 

ONE-PHASE 

a G 

TWO-PHASE 

a o,G 

80 

894 

3.75 

177 

163 

100 

613 

16.8 

183 

167 

120 

338 

49 

176 

139 


The discontinuity is consistently weaker on the vapor side, with the 
strongest discontinuity occurring on the saturated liquid side at a 
temperature of 80 K where a Q L is an incredibly small 3.75 meters per 
second. Normally, the boundary sound speed associated with the single- 
phase side is the one recognized as the actual boundary acoustic velocity. 
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For example, when using the saturation tables, in reference 6, it is tacitly- 
understood that these boundary values are to be taken as the normal single- 
phase values. In the case of water, this discontinuity is discussed in 
reference 7» and a Q is plotted isothermally as a function of void fraction 
in figure 1 of reference 8 and again shows liquid phase values close to 
zero . 

In addition to the (IHE) model, various non-equilibrium models have 
also been developed. One of these introduces a non-equilibrium feature 
by neglecting changes in phase composition during passage of the sound 
wave. As with a^, the basic sound speed expression used is 


a = [-v 2 /(3v/9P) 

s 


1/2 

] 


( 1 ) 


where 


v 


(l— x) v L + X V Q 


( 2 ) 


with x being the vapor mass fraction (quality) and v T and Vj, the 
liquid and vapor specific volumes. While both the (IHE) and non-equilibrium models 
assume constant mixture entropy for the wave propagation, the latter model also 
assumes constant liquid phase entropy as well as constant gas phase entropy; that 
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is, s, s L , and s G are all individually constant during the passage of the 
sound wave as is the quality. With these assumptions equation (2) is differ- 
entiated into 

(§) s = Cl ' x> (l?) s>L + x (t?) o>(! (3) 

The partial derivatives on the right hand side of this equation are 

2 2 2 2 

respectively replaced by (- v T /a T ) and (- v n /a n ), and the resulting 

L L It vjr 

expression is substituted into equation (l) to give the metastable, 
triply isentropic, (M, 3s) sound speed as 


a 3 = v/{(l - x) ( VjV a ^) 2 + x(v G /a G ) 2 } 


(b) 


Here, a-^ and a^ are the isentropic, single -phase boundary values. This 
equation is used by Karplus (reference 7)5 Wood (reference 8), and Dvornichenko 
(reference 9) in a form that substitutes vapor void fraction for quality. 

In reference 10, Wallis generalizes the isentropic phase constraints used 
in defining a^ and a^ to permit other acoustical vibrational constraints 
such as isothermal phase propagation to be introduced into the theory of sound 
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propagation. When these constraints are the same as in this section, 
his equation (.2.1+6) is identical with equation (4) above, and his 
equation (2.50), with the void fraction replacing the quality, is 
identical with the equations derived in references 7~9* 

The important features of this (M,3s) model are summarized below; 
(l) In contrast to the (IHE) model, the acoustic velocity of the 
(M,3s) model reduces to the usually measured single-phase velocities 
for x = 0 and x = 1. With a Q used to denote the sound speed 
associated with the (IHE) model, and with a^ used to denote (M,3s) 
sound speed, the boundary value comparison in meters per second for the 
two models using nitrogen at 100K gives 


a„ T = l6 . 8 


0 ,L 

3,L 

a 0,G = 167 

a 3,G 


613 

183 


(2) When sound speed is plotted against quality at constant 
temperature and pressure using the (M,3s) model, the resulting curve 
exhibits a minimum with the minimum sound speed being smaller than either 
boundary value. For nitrogen at 100 K the minimum velocity of sound is 
about 77 meters per second and occurs very close to the saturated liquid 
boundary at a quality of about 0.049. (As will be shown later, the (IHE) 
model predicts a monotonic curve.) Measured data for the two-component 
air-water mixture also shows a similar minimum (reference 7 and reference 
10, p. 266) and a similar pattern holds for measurements on a single- 
component steam-water system as shown in reference 11. 
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While qualitatively correct, the (M,3s) model overpredicts the mini~ 
mum sound speed in the air<-vater system of Wallis (reference 10, figure 9.12) 
by approximately 20% . In order to account for this discrepancy, Karplus 
(reference 7 ) first transforms equation ( 4 ) to replace quality by void 
fraction and then makes an ad hoc modification in the (M,3s) model by 
replacing isentropic propagation in the phases by isothermal propagation; 
that is, he permits heat transfer between the phases. Actually Karplus 
did not concern himself with the liquid phase isothermal modification 
because two further approximations, p L » p Q and p L » p Q a^, enable him to 
transform equation (1|) into 



( 5 ) 


where a is the void fraction. In this form the question of 

isentropic versus isothermal propagation for the individual 

phases is not yet addressed; however, since all reference to 

the liquid phase velocity of sound has been eliminated, no 

isothermal substitution for isentropic propagation need be made for 

the liquid. For atmospheric air at room temperature an ideal gas 

substitution can be made for the gas phase speed of sound according to 
2 

a Q = yRT if isentropic, or simply RT, if isothermal, Choice of the latter 
then yields the following equation after using P = pRT, 
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>1 


( 6 ) 


2 

a g x P[p L a(l - a)] 


for an isothermally modified model (M,s,2T) for which a^ now represents 
isothermal sound propagation in both phases even though the approximate 
equation (6) only involves the gas phase sound speed. The single s 
is retained in the classification notation on the assumption that 
Karplus models the mixture as isentropic. Without using the above two 
Karplus approximations concerning the magnitudes of sound speeds and 
densities, the (M,s,2 T) sound speed may also be formed from equation 
(4) by the appropriate isothermal substitutions of a T m for a T 

Lj , 1 -Li 

2 

and a^ ^ for a^. Then, with these new symbols representing {-v (3P/3v)^} 
respectively evaluated at the saturated liquid and vapor boundaries, 

a 2 = t/{(1 - x)( V a L,T )S + x( V a 0,T )2! <7) 


The significant results of this modified, ad hoc, (M,s,2T) theory- 

are: 

(1) It correctly predicts the magnitude of the sound speed in the minimum 
point region in the Karplus experiments (reference 7)» which are also shown in 
Wallis (reference 10, p. 2 66 ). 

(2) Equation (7) for the isothermal, two-phase, sound speed does 
not yield the usual single-phase, isentropic sound speeds at the boundary 


11 



points x - 0 and x = 1, or a - Q and a - 1. (The approximate equation (6) 
gives boundary- velocities that are both infinite, but Karplus indicates 
negligible error in the use of equation (6) for 0.002 < <x < 0.94.) 

The difference between isentropic and isothermal propagation is 
negligible for pure-phase water at room temperature because the ratio 
of heat capacities is close to unity. However, pure-phase air, treated, 
as an ideal gas with a heat capacity ratio of y = 1.4, will have a ratio 
of isentropic to isothermal sound speed of /l. 4 = 1.183 and, coincidentally, this 
is about the same order of magnitude difference as the original overprediction of 
the minimum sound speed inherent in the (M,3s) model. The features of the 
air-water system are shown in figure 1 where the sound speeds a^ and 
a^, corresponding respectively to the (M,3s) and (M,s,2T) models of 
equations (4) and (7) are plotted. Discussion of the third curve, a^, 
will come later. In the computation, air was treated as an ideal gas and the 
water data were taken from reference 12. Wallis in reference 10, page 266, 
shows experimental data for the air-water system at a similar temperature along 
with the corresponding theoretical lines. A reproduction of his figure 9.12 is 
shown herein as figure 2. The line labeled by Wallis as "isothermal" corresponds 
to the (M,s ,2T) model, ag, and his line labeled "adiabatic" corresponds to the 
(M,3s) model, a^. 

Because of the steeply rising curves on the gas side of figure 1 
where the void fraction is close to one, the difference in the gas 
boundary values of the sound speeds for isothermal and isentropic propa- 
gation is obscured. In order to illustrate these differences better and 

also to emphasize the mass dominance of the liquid phase in the vicinity 
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of the minimum point, the void fraction, a, is transformed to the quality, x, 
by use of 

x - av^/{vg a(v & - v L )} (8) 


Then the same sound speeds, a^ and a^, are replotted in figure 3. 

Now the x = 1 values of a^ and a^ are clearly evident, as is the 
mass dominance of the liquid phase. The minima in figure 1 show a rough 
balance between liquid and vapor volumes $ that is, the value of the void 
fraction that minimizes any one of the curves is about 1/2. Indeed, 
the approximate form of a 2 from equation (6) has its minimum exactly 
at a void fraction of l/2. In figure 3, the previous, apparent balance 
between the phases is completely eliminated as the minima are crowded 
against the pure liquid axis. This reflects the enormous density 
difference between water and air. Again, discussion of the a^ curve 
is deferred. 

Further complications in matching theory and data are introduced 
by Semenov and Kosterin in reference 11, In contrast to the data shown 
by Wallis, their measurements on single-component steam-water systems 
give results in apparent quantitative agreement with the basic, unmodified 
(M,3s) theory. At first, their results may suggest a fundamental dis- 
tinction between single- and two-component two-phase systems. On the 


13 



other hand, the uncertainty in their acoustic velocity measurements was 
reported to be +9%> and this may be sufficiently large to obscure the possible 
differences between the data and the (M,3s) theory* 

In view of these conflicting experimental results and theoretical models, 
a new metastable theory of sound propagation in two-phase systems was sought 
which would predict the correct minimum in the acoustic velocity-void fraction 
curve as does the (M,s,2T) model and would achieve isentropic single-phase 
velocities at x = 0 and x = 1 as does the (M,3s) theory. The new approach 
retains two key concepts of the (M,3s) theory, - that a local acoustic 
disturbance is again assumed to propagate in a homogeneous two-phase mixture, 
(l) without phase change, 

(.2) at constant mixture entropy. 

These postulates make the metastability similar to isentropic expansion of a 
super-heated vapor in a nozzle to a pressure less than saturation without 
condensation occuring immediately. In this well known case, only the gas 
phase is present at the onset of the metastability; whereas, metastable 
propagation of sound in a homogeneous two-phase mixture has the added 
complexity of requiring a similar non-equilibrium effect for each phase at the 
onset of metastability. At this point the new theory drops the (M,3s) 
additional requirements of constant s^ and s^ for a different concept of 
metastability described in the next section. 

METASTABILITY 

In order to develop a new description of two-phase metastability 
from the stated postulates, the total differential of the 
entropy function for a single phase is taken from 



thermodynamics and written for each phase as 


dS L 



dT - (fj) dP 
di P,L 


(9) 


and 


ds G F dT - < tF ) P;0 ap 


( 10 ) 


For a saturated two-phase system of one component , these two 
differentials are expanded about two different states, each at the same 
temperature and pressure. One of these states corresponds to the satu- 
rated liquid, and the other, to saturated vapor. Let these states be 
the ones shown as points L and G in the sketch of figure Ma) which 
shows a pair of isobars at P and P+dP passing from the liquid region 
into the two-phase region and finally into the vapor region. L' and G' 
are respectively saturated liquid and vapor states at the lower isobar. 
Starting from the point L and assuming an isentropic expansion occurs 
which drops the pressure (dP < 0), total phase equilibrium - including 
phase change- requires the expansion to terminate at the two-phase point 



marked as A. Similarly, an isentropic acoustic expansion from G would 
terminate as the two-phase state B. 

In order to impose the condition dx = 0 during the wave passage, 
penetration into the two-phase region must be avoided. To make this 
possible, a second sheet of me tas table -unstable properties is assumed 
to underlay the normal two-phase equilibrium property sheet. For pur- 
poses of discussion it is sketched in figure 4(b) in a manner similar 
to a van der Waals metastable-unstable region. Here the usual equilibrium 
isobars are shown as solid lines P and P+dP with L, L* , G, G' being the 
same former saturation states. The dashed curves represent one-phase, 
metastable-unstable continuations of these isobars. Now the passage of 
a disturbance (dP < 0) from the isothermal and isobaric points L and G 
without phase change shifts these points to the new isobar P + dP such 
that the final phase points have the same (and new) final temperature. 

For the way figure 4(b) is sketched it is impossible for the individual 
phase points to shift about local isentropes passing through L and G and 
satisfy the final isothermal condition. If this sketch were a true 
representation of thermodynamic properties, the (M, 3s) theory of sound 
propagation could not possibly be valid. On the other hand, it is possible 
of course to alter the sketching of the isobar P + dP such that the pair of 
local isentropic expansions emanating from L and G end on the P + dP isobar 
with the same temperature in conformity with at least this portion of the 
(M,3s) theoretical requirements. Obviously, the thermodynamic pattern of 
properties is rigidly set, and any theory of two-phase sound propagation 
must conform to this pattern. Later, such a conformance test will be 
derived for the (M, 3s) model which, when tested for several substances 
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including nitrogen, indicates that this model is invalid with respect to 
the thermodynamic data used. For the present, it suffices that the two 
sketches indicate the potential metastable thermodynamic states which can be 
associated with two-phase sound propagation. 

DERIVATION OF THE NEW EQUATION 

In the derivation of the new equation for metastable sound propaga- 
tion the basic starting point is again equation (l), and this means finding 

a way of evaluating (3v/3P) .. Likewise, the mixture specific volume continues 

s 

to be described by equation (2), and it is perturbed by the passage of the 
wave to give 


dv = (l - x) dv 

Li 


+ x dv G 


( 11 ) 


since propagation is again assumed to occur without phase change. The 
single-phase specific volumes are ordinarily functions of temperature 
and pressure whose total differentials are 


dv L = 'If' 


dT + 


P,L 


(— ) 


dP 


T,L 


( 12 ) 
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and 



(— ) 
V 9T ; 


dT + 


P,G 


(— ) 

'3P' 


dP 


T,G 


(13) 


where dT and dP are the same in both equations and are the perturba- 
tions associated with the sound wave. Substitution of the last two equations 
into the preceding one gives 


dv = ( It^ p dT + ( iP T dP 


(lb) 


where the mixture quantities (3v/3T) p and (3v/3P) T are given by 


# p ■ & - 


X) <fe) - *(S) 

3t p,l 9t p,g 


(15) 


and 


(|f) =(l-x)(& ■ 

3P T 3P T,L 


/3vs 

x( 3P> 


T,G 


(16) 
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For later use, the third partial derivative of the equation of state is 
formed in the usual way to give 


( 2£) = „ (it) / (li) 

V 3T ; p ' ^9P ; t 


(IT) 


Since the problem is to form (3v/9P) * the mixture entropy is now 

s 

introduced as 


s = (l - x) s L + X Sq 


(18) 


and differentiated without phase change to give 


ds 


(1 


x) ds L 


+ X dS G 


(19) 


The quantities ds^ and ds^ are eliminated by use of equations (9) and (10) 
to give 





(20) 
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( 21 ) 


where Cp is defined for the mixture as 



X> °P,L + X °P,G 


and (3v/3T) p has been previously defined in equation (15 ). It is noted 

that the heat capacity is the usual T(3s/ 3T) p and is not infinite in 

the metastable theory. This is contrary to its infinite behavior in 

stable two-phase theory for pure substances. 

Now all that is required to form (3v/3P) is to change the 

s 

independent variables in equation (20 ) from dT and dP to dv and dP. This 
is accomplished by rewriting the equation of state by substituting 
equation (17) into equation (l^). This gives 


dT 


(— ) 
'3P' 


dP + 


M) 

vv . 


dv 


( 22 ) 


and this is used to eliminate dT in equation (20) with the result that 


ds = 


(3T) 

l 3v'. 


dv + 




m . 

*3p' 

V 


(— ) 

V3T ; . 




dP 


(23) 
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From this, the partial derivative (3P/3v) 

s 


is read as 



C (— ) 

°P 3v P 

n /3T» m /3V\ 

C P 3P ~ T 3T p 
v P 


(2k) 


Multiplying numerator and denominator by (3P/3T)_^ and then using equation 
(17) in the numerator gives 



c A 

C P ^3v' 


T 


C - T(— ) (— ) 

°P M 3t' '3T ; _ 
v P 


(25) 


In stable one-phase thermodynamics the denominator would be the heat 
capacity at constant volume. To test this, the equation of state is 
rewritten as 


dP = 


( 9P) 

V 3T ; 


dT + 


M) 

( dv\ 


dv 


(26) 


and this is used to eliminate dP in equation (20). The result is 
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( 27 ) 



(3V) (9P) 

v 3t' V 3T ; 

P v. 


.} dT - ( £> T 


( i> p ^ 


From its basic definition, 


C v ■ 1 <H> 


(28) 


Therefore , it follows from the last two equations that 


C = C — T (— — ) (-5E-) 

v L P T 9T „ W 

P v 


( 29 ) 


and, from equations (25) and (29) 


( 3P) 


/3P-, 

■ y <37> 


T 


( 30 ) 


where 


Y = C p /C 


v 


( 31 ) 
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When equation (30) is substituted into equation (l), the final form of the 
acoustic velocity, now designated as a^, becomes 


? 1/2 

a 1 = [ -Y v /(3v/3P) T l 


(32) 


In summary, the expression for a^ in equation (32) represents the new 
metastable equation. It has assumed no phase transition during an acoustic 
disturbance and has assumed that the total entropy of the system (not of the 

individual phases) remains constant. As will be shown, the values of a^ 
reduce to the respective single-phase values of sound speeds in the limits 
of all liquid or all gas. 

For computational purposes one must know the ordinary single-phase values 


° f C P,L’ C P,G’ C V,L’ C V,G’ ( 9v / 9T )p >L > ( 9v / 9T )p }G > (3v/3P) T ^ L , ( 9 v / 9P )t,G’ 
v T , and in order to obtain a sound speed versus either quality or void 

fraction curve at a given temperature and pressure. In conjunction with 
this, it is noted that the right side of equation (2l) reduces to the 
ordinary single-phase heat capacities at x = 0 and x = 1. A similar 
reduction occurs at these qualities for (3v/3T)p and (3v/3P)^ as is 
evident in equations (15 ) and (l6). Substitution of equation (IT ) in 
equation (29) gives an alternate form of the constant volume heat capacity 
as 


C v = C p + T(3v/3Ty (3 v/3P) t 


(33) 
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At x = 0 this reduces to the usual single-phase form 


C = C + 
V,L P,L 


T(3v/3T ) / (3v/3P) 

iT 5 i-J X 5 i-J 


(34) 


A similar reduction occurs at x = 1, Finally, equation ( 8 ) or the 
alternate form 

a = x v^/v (35) 

can he used to convert quality into void fraction or vice versa. 

RESULTS AND DISCUSSION 
Two-Component System 

Referring hack to figures 1 and 3 for the air-water mixture, the 
remaining curve to he discussed is the one labeled a^ and is, of 
course, from the new theory embodied in equation (32). The first figure 
shows that it is practically coincident with the ad hoe (M,s,2T) curve, a^, 
which, according to the results in Wallis (ref. 10) correlates the 
available experimental data about the minimum point, In figure 3 the 
a^ curve subsequently deviates from the a^ curve well past the minimum 
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point and intersects the curve at the pure gas side as both the a^ and 
a^ curves satisfy the gas boundary condition of isentropic propagation. 

An insight into the mechanism explaining similar behavior for a^ and ag 
at the minimum point and liquid boundary and different behavior at the gas 
boundary is obtained from figure 5, which plots the two-phase heat capacity 
ratio versus the void fraction. Bearing in mind that the minimum in the 
previous plots is located approximately at a = £-, this figure shows how 
the two-phase heat capacity ratio adjusts itself to an approximate value of 
unity in the minimum region from a void fraction of zero to about 0.95 
before rising steeply to the pure gas value of y = 1.1+ at a = 1. The approxi 
mate value of unity is needed for correspondence with isothermal sound 
propagation postulated in the (M,s,2T) model to fit the experimental points, 
and the 1.4 value is necessary to satisfy isentropic propagation in the pure 
gas phase. 

One-Component Systems 

With the new equation accomplishing the dual objective of correlating 
the minimum and boundary points for the air-water system, a comparative 
study was made for parahydrogen and nitrogen, for which similar experi- 
mental data Were unavailable. Both substances are in the cryogenic domain 
with nitrogen of primary interest because of its future role in the NTF 
program, and both are single-component two-phase systems. While steam 
and water also form a one-component system, this two-phase mixture will not 
be included in the comparative study. Instead, it is used to develop a 
thermodynamic criterion for determining if the (M, 3s) theory can be correct 
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as implied by Semenov and Kosterin in reference 11 for steam and water* 
Nitrogen 

In the nitrogen study three plots similar to those shown for the air- 
water system were constructed for 80K, 100K, and 120K. These temperatures 
were arbitrarily selected to represent nitrogen’s entire two -phase 
temperature range, 63.1^8 K <_ T <_ 126. 2K. In addition, the (IHE) model 
was used to show a fourth velocity of sound designated as a Q . Cryogenic 
tunnels such as NTF will normally have a^l and x >0.95* 

At the nitrogen temperature of 80K, which is furthest from the 

critical point, figures 6, 7> and- 8 most closely resemble the three 

corresponding figures for the air-water system. The chief difference between 

one- and two-component systems as far as these figures are concerned is 

the y value. For the air-water system it is close to unity while for 
Li 

nitrogen it is slightly greater than 2.1. This is shown in figures 5 
and 8. The latter figure also shows the very rapid adjustment that the 
two-phase y makes to values very close to unity for a large range of 
a (although not as large a range as that for the air-water system). 

Because of the rapid adjustment in y, the minimum region of figure 6 
shows a.^ and a 2 as being essentially one curve. 

In figure 6 the (IHE) velocity of sound previously defined as a Q 
is shown to be smaller than the minimum value of a^ for a very large 
range of a before intersecting a in the vicinity of a =1. This 
region is difficult to visualize because of the convergence of four 
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steeply rising curves; therefore, these curves are transformed by using 
the quality as the independent variable in figure 7 • Here it is seen 
the a Q curve intersects a^ at an approximate quality of 0.3, and 
subsequently makes an apparently tangential contact with the a^ curve at 
a quality of approximately 0.55* This tangential contact raises an 
interesting question about the meaning of points which appear to be 
simultaneously stable and metastable . 

When the temperature of the nitrogen is changed to 100 K, the pattern 
associated with the previous three figures begins to shift a little 
although the basic features remain the same as seen in figures 9, 10, and 
11. The a^ and a^ curves merge in a narrower minimum region of figure 9 
with the minimum sound speed still being less than that predicted by the 
(M,3s) model. The narrowing of the merged region is accompanied by a more 
visible shift in the y pattern as indicated by comparing figures 8 and 11. 
No longer is y equal to one over a large fraction of the void fraction 
axis as in figure 8, even though it is still very close to one in its flat 
minimum region. Finally, the prior interaction and tangency characteristics 
of the equilibrium sound speed, a Q , persist. 

The variation in the triple pattern of plots becomes stronger when the 
temperature is shifted closer to the critical value for the final value 
of 120 K as seen in figures 12, 13, and l4. Now there is no region in figure 
12 where the a^ and a^ curves merge. This strong change is accompanied 
by an equally strong shift in the y pattern of figure 14 where no region 
exists for which y can be approximated by unity. Therefore the (M,s,2T) 
model should not accurately predict the minimum region in this case 
because the necessary assumption of isothermal propagation is no longer valid. 


27 



Finally, the intersection and tangency characteristics in the a^ curve 
persist, as is evident both in figures 12 and 13. 

Parahydrogen 

Paraliydrogen has triple and critical temperatures of 13.803 K and 

32.976 K (reference 13). From this range only the single temperature of 

25 K was selected for the same kind of plots as before. This gives a 

T 

reduced temperature of zr~ = 0.758, which is very close to the reduced 

J-c 

temperature of 0.792 for nitrogen at 100 K. From a principle of 
corresponding states view, it would be anticipated that figures 15, l6, 
and 17 for parahydrogen would.be very similar to those for nitrogen as 
shown in figures 9» 10, and 11. This turns out to be only partially 
true. Figure 15 shows that while the a , and a g curves are very close together 
there is no region of overlap as is evident in the corresponding nitrogen 
figure 9- Subsequent comparison of figures 17 and 11 shows that y for 
the parahydrogen has a narrower flat region about its minimum point which 
leads to a poorer approximation to the y = 1 condition needed to merge the 
a^ and ag curves in the vicinity of their minimum points. 

Steam-Water 

As a concluding part of this investigation, attention is now turned 
to the steam-water experiments of reference 11 which apparently support 
the triply isentropic model (M, 3s) within the context of an uncertainty 
in sound speed measurement of +9%. The triple plots used in the compara- 
tive study of nitrogen and parahydrogen will be replaced by an analytical 
treatment involving the use of equations ( 9 ) and (10). 
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Since pressure and temperature equilibrium is assumed in the 
metastable models, the perturbation differentials dP and dT are the same 
for the liquid phase, the vapor phase, and the two-phase element. In the 
(M, 3s) model ds = ds^ = ds^ = 0. Applying these assumptions to equation 
(9) and (10) gives the perturbation pressure-temperature ratio dP/dT 


dP _ ^P,L ,3T . 
dT " T l 3v' 

r 


(36) 


and 


dP _ C P, G/3T, 
dT' T 


(37) 


Equating the two expressions for dP/dT gives 


C M) 
C P,L l 3v ; 


= C. 


(S 


P,L ^ G ' 9v P,G 


( 38 ) 


which will be used in a nondimensional ratio form, r, by defining 
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r = 


(39) 


P /9T\ /3ts 
C P,G 9t P,L 9T P,G 


where, of course, 


r = 1 


(40) 


This result is important in that the postulates of the (M,3s) 
model predict a property law which can be checked by using suitable 
tables of thermodynamic properties. For the steam-water system it was 
checked by using the properties found in table 1 of reference 11, and by 
transforming r to accommodate them. This is accomplished by rearranging 
the thermodynamic single-phase equation for the difference in specific 
heat capacities to 


(2z) = \ / c g c v 77T^ Y 
9t P y ■ • * i ▼ 


(41) 


Applying this successively to the liquid phase and to the gas phase and 
then substituting both into equation (39) gives 
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Introduction of the tabulated data from reference 11 gives the following 
results for steam-water at four saturation temperatures: 


T(K) 

r 

452.2 

746 

470.5 

461 

484.5 

318 

522.3 

112 


Since r fails to satisfy equation (40) the conclusion, based on 
the equations developed in this paper,, is that the (M, 3s) theory is 
invalid and that the +p% uncertainty in the steam-water sound speed mea- 
surements of reference 11 is sufficient to conceal this. 

In the case of the air-water system it is known that the (M, 3s) 
model does not fit the experimental data. From the above analysis it 
would be expected that r 4 1 and indeed the computation gives r = 84300. 

For parahydrogen experimental data is lacking for the required sound 
speed curve. Falling back on the r computation gives r = 39-3 which 
invalidates the (M,3s) model in the context of the theory of this paper. 
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A similar conclusion is drawn for nitrogen where, again, experimental 
data for the sound speed curve are lacking; for, at the three temperatures 
80K, 100K, and 120K, the respective values of r were calculated as 58 O , 

64.4, and 7*80. 


SUMMARY OF RESULTS 

A new method of calculating speed of sound for two-phase flow has 
been presented. The new equation assumes no phase change during the 
propagation of an acoustic disturbance and assumes that only the total entropy 
of the mixture remains constant during the process. It is shown that the 
new equation predicts single-phase values for the speed of sound in the 
limit of all gas or all liquid and agrees with available two-phase, air- 
water sound speed data. Other expressions used in the two- 
phase flow literature for calculating two-phase, metastable 
sound speed are reviewed and discussed. Comparisons are 
made between the new expression and several of the previous 
expressions — most notably a triply isentropic equation as used, among 
others, by Karplus and by Wallis. Appropriate differences are pointed out 
and a thermodynamic criterion is derived which must be satisfied in order 
for the triply isentropic expression to be thermodynamically consistent. 

This criterion is not satisfied for the cases examined, which included two- 
phase nitrogen, air-water, two-phase parahydrogen , and steam-water. Con- 
sequently, the new equation derived is found to be superior to the other 
equations reviewed. 
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The importance of this work to cryogenic wind tunnels is that is 
properly defines the two-phase speed of sound to use as the mixture approaches 
the limit of all gas. In particular, the correct speed of sound approaches 
the single-phase gas velocity of sound and not a value substantially 
less as predicted by an isothermal formulation. In other words, in regions 
of small amounts of liquid mixed with gas (vapor mass fraction greater than 
0.98)j the speed of sound will be very close to that value expected for 
single-phase flow. 
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Air -water 



Figure 1.- Different sound speeds for an air -water system at 290 K as 
a function of void fraction. 
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Figure 3.- Different sound speeds for an air -water system at 290 K as 
a function of quality. 
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(a) Equilibrium description. 



(b) Metastable description. 

Figure 4.- Description of metastability on thermodynamic properties. 
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Figure 5.- Ratio of specific heats for an air -water system at 290 K 
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Figure 6 



Different sound speeds for two-phase nitrogen at 80 K as 
a function of void fraction. 



Figure 7.- 



Figure 8.- Ratio of specific heats for two-phase nitrogen at T = 80 K. 




Figure 9.- Different sound speeds for two-phase nitrogen at T = 100 K as 
a function of void fraction. 



Figure 10.- Different sound speeds for two -phase nitrogen at 100 K as 
a function of quality. 








Figure 14.- Ratio of specific heats for two-phase nitrogen at T * 120 K. 





Figure 16.- Different sound speeds for two -phase parahydrogen at T 
as a function of quality. 
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r Figure 17.- Ratio of specific heats for two -phase parahydrogen at T 
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= 25 K. 
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